Ab initio properties of Li-Group-II molecules for ultracold matter studies 
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We perform a systematic investigation of the electronic properties of the 2 E + ground state of Li- 
alkaline-earth dimers. These molecules are proposed as possible candidates for quantum simulation 
of lattice-spin models. We apply powerful quantum chemistry coupled-cluster method and large 
basis sets to calculate potential energies and permanent dipole moments for the LiBe, LiMg, LiCa, 
LiSr, and LiYb molecules. Agreement of calculated molecular constants with existing experimental 
data is better than or equal to 8%. Our results reveal a surprising irregularity in the dissociation 
energy and bond length with an increase in the reduced mass of the molecule. At the same time 
the permanent dipole moment at the equilibrium separation has the smallest value between 0.01 
a.u. and 0.1 a.u. for the heaviest (LiSr and LiYb) molecules and increases to 1.4 a.u. for the 
lightest (LiBe), where 1 a.u. is one atomic unit of dipole moment. We consider our study of the 
2 E + molecules a first step towards a comprehensive analysis of their interactions in an optical trap. 



I. INTRODUCTION 

There is growing interest in molecules formed by one 
alkali-metal atom and one alkaline-earth or rare-earth 
atom [ItiTT]. In contrast to group I or group II di- 
atomic molecules, these molecules have an unpaired elec- 
tron that allows for manipulation of the X 2 S + ground- 
state molecule with both electric and magnetic exter- 
nal fields. Moreover, it has been suggested that 2 E + 
molecules are good candidates in which to explore con- 
trolled chemical reactions at ultracold temperatures [H] . 
The unpaired spin in the ground state of these molecules 
should provide a new handle for the control of the reac- 
tion dynamics based on the interplay between intramolec- 
ular spin-rotation couplings. In addition, when placed 
in an optical lattice, these polar molecules can inter- 
act with each other via both electric dipole-dipole and 
magnetic spin-spin forces. It is therefore possible to en- 
gineer unusual forms of interactions for quantum sim- 
ulation of lattice-spin models and topological quantum 
computing Q] , controlled preparation of many-body en- 
tangled states [1], and high-precision measurements of 
fundamental constants [13] . 

Another interest in Li-group II molecules stems from 
prospects to achieve optical Feshbach tuning of scattering 
properties in ultracold gases of the individual atoms with- 
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out substantial loss of atoms from a trap. This photoasso- 
ciative tuning, pioneered for homonuclear molecules |14L - 
121] . becomes possible due to the existence of long-lived 
excited molecular states near narrow intercombination 
lines of the alkaline-earth or rare-earth atoms. In turn, 
this might enable efficient ways to form gases of polar 
molecules without substantial loss by photoassociation. 
In essence, a two-photon optical Feshbach resonance can 
be used to couple two colliding atoms to a vibrational 
level of the molecular ground state. The suppressed ef- 
fect of excited-state spontaneous decay makes efficient 
coherent molecular formation possible. Knowledge of the 
electronic and ro- vibrational properties of the molecules 
will help to find optimal pathways. Characteristics such 
as the permanent dipole moments of the ground state 
will determine the anisotropic interactions between these 
molecules in an optical trap. 

Some key theoretical predictions of the electronic 
structure and dipole moment of the ground and low-lying 
excited states of LiYb [9] and for the ground state of 
LiSr [5] have recently been obtained. Other work has 
studied RbYb [2] and RbSr |5]. However, the physi- 
cal origin of the chemical bond of these mixed species 
is not yet fully understood. Continued theoretical ad- 
vances are necessary to unravel the properties of these 
molecular systems. To the best of our knowledge, only 
the LiMg and LiCa ground states were experimentally 
investigated [HI E21 EH] , which allows us to compare our 
predictions with experimentally obtained molecular pa- 
rameters |lll 123) . There are several experimental groups 
working to achieve the formation of the LiYb molecules 
at ultracold temperatures [3 [10] . An understanding and 
quantitative description of these molecules might help to 
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define conditions under which these ultracold molecular 
systems can be created. 

Our theoretical study is devoted to a systematic in- 
vestigation of the X 2 £ + ground state properties of the 
polar Li-alkaline-earth dimers such as LiBe, LiMg, LiCa, 
LiSr, and the rare-earth dimer LiYb. These are inter- 
esting and challenging systems as the valence electrons 
of the atoms form a relatively weak covalent bond [2"5] . 
The alkaline-earth atom has a closed outer electron shell 
and a Hartree-Fock interaction energy results in a purely 
repulsive potential curve. We use the unrestricted (U) or 
partially spin-restricted (R) coupled cluster (CC) method 
with either single and double excitations augmented with 
pertubatively calculated triples (CCSD(T)) or with sin- 
gle, double, and triple excitations (CCSDT) to treat the 
correlations of core and valence electrons. All approaches 
are used in conjunction with correlation-consistent basis 
sets. The dipole moment is calculated using a finite field 
approach with a four-point numerical differentiation for- 
mula. The calculations are performed with CFOUR [26] 
and MOLPRO [37] suites of programs. 

Presumably, the CCSDT method is more accurate in 
the description of the potential energy surface than the 
CCSD(T) method, where triple excitations are treated 
perturbatively. The results of the unrestricted (U) 
method are more accurate than that of a restricted 
(R) calculation. The CC method is based on a single- 
reference configuration and is not suitable for molecules 
where multiple-configurations are relevant. A practical 
indicator of how well the CC method performs is based 
on Tl diagnostics [28]. Our calculation shows that Tl 
is small and varies from 0.015 to 0.035. This indicates 
that a single-reference method is valid for these dimers. 
In fact, a small Tl for all interatomic separations implies 
that the long-range interaction potentials as well as the 
dissociation energies can be accurately described. 

Our results include the ground-state potential curve, 
bond distance, harmonic frequency, dissociation energy, 
and permanent dipole moment for LiBe, LiMg, LiCa, 
LiSr, and LiYb. For completeness, we calculate vibra- 
tional energies of the ground-state potentials. Due to 
limitations in computational resources we were not able 
to obtain results for all methods, basis sets, and Li-X di- 
atoms. The most accurate potential energy curves and 
dipole moments are presented in the figures and tables 
with the exception of Tab. |Vl] which compares molecular 
constants obtained with the various methods and basis 
sets. 

All potentials are corrected for basis set superposition 
errors (BSSE) using the counterpoise procedure of Boys 
and Bernardi |33| . The sign of the dipole moment de- 
pends on the position of the origin. We apply the fol- 
lowing convention throughout the paper: The alkaline- 
carth/rare-earth atom is positioned at the origin and the 
Li atom moves along the positive x-axis. A positive sign 
of the dipole moment corresponds to the charge transfer 
from the Li atom towards the alkaline-earth/rare-earth 
atom. 



A. LiBe molecule 

There exist several theoretical calculations of the 
X 2 E+, 2 electronic ground state of the LiBe molecule. 
Jones j29| used a density functional formalism based on 
a local density approximation to predict molecular con- 
stants of the LiBe, LiMg, and LiCa molecules. Their 
dissociation energy is D e /(hc) = 3705.9 cm -1 , while the 
equilibrium bond length R e is 4.9cto and the vibrational 
constant uj e /(hc) = 340 cm -1 . Such strong bonding of 
the LiBe molecule was explained by substantial electron- 
charge transfer between s and p orbitals. More advanced 
MRCI calculations of the low-lying doublet and quar- 
tet states of LiBe was performed by Fischer et al. [30j . 
In this calculation the ground state is bound by 2497.5 
cm -1 , R e = 4.95eto, and Lu e /(hc) — 302 cm -1 . The per- 
manent dipole moment is d e = 0.95 a.u. at R e . Another 
MRCI calculation [35] using all single and double excita- 
tions obtained D e /(hc) = 2336.35 cm -1 , R e = 4.927ao, 
and ui e /(hc) = 300 cm" 1 , which are close to those of 
Ref. [30]. However, the dipole moment d e — 1.33 a.u. 
of Ref. [25] differs significantly from that of Ref. [30] , 
The effect of inner-shell correlations on the molecular 
electronic structure was investigated. Later an ab initio 
configuration-interaction study was reported by Marino 
and Ermler [3T]. Their value of the dissociation energy 
D e / (he)— 2014 cm -1 is smaller than that in previous 
studies. In summary, the previous theoretical calcula- 
tions disagree on dissociation energy and dipole moment 
values but have much better agreement on the value of 
the bond length. Our goal is to test results of the previ- 
ous predictions by performing a coupled cluster calcula- 
tion accounting for high-order correlation effects. 

In our calculation the UCCSD(T) method is used 
to obtain the ground state potential and dipole mo- 
ment of LiBe as a function of interatomic separa- 
tion R. These characteristics are calculated with 
the large basis set (aug-cc-pV5Z-DK) of Prascher et 
al [35], which includes scalar relativistic effects ac- 
counted for by the Douglass-Kroll Hamiltonian. Atomic 
basis sets include (15s8p5d4f3g2h) primitive Gaussian 
and [7s6p5d4f3g2h] constructed functions for Li and 
(15s9p5d4f3g2h)/[7s6p5d4f3g2h] for Be. 

Results of our calculation are shown in Figure [l] Com- 
paring our data with results available from other theo- 
ries and experiment we conclude that R e — 4.87ao is 
close to the experimental measurement of R e — 4.894ao 
[52] . while the well depth D e /(hc) = 2406.89 cm" 1 is in 
best agreement with the prediction of [25]. In Ref. [32] 
no value of D e is given. The dipole moment value d e 
— 1.41 a.u. at the equilibrium separation is even larger 
than that of [25] , The permanent dipole moment aver- 
aged over the ground state vibrational wave function is 
equal to 1.37 a.u. It is slightly lower than the value of the 
dipole moment function at R — R e . The LiBe molecule 
has the largest dipole moment of the molecular systems 
considered in this work. 

Bound states of 7 Li 9 Be are calculated using a discrete 
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FIG. 1: Top panel: electronic 2 E + ground-state potential of 
LiBe; bottom panel: permanent dipole moment function of 
the ground state of LiBe. The vertical dashed line indicates 
the equilibrium separation of the ground state potential, at 
which the dipole moment equals 1.41 a.u. The vibrationally- 
averaged dipole moment of the v = level is equal to a slightly 
smaller value of 1.372 a.u. 



TABLE I: All vibrational energies of the 7 Li 9 Be X 2 E+ ground 
state potential for the rotational state N = 0. 
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Energy/(/ic) (cm 1 ) 




Energy/(/ic) (cm 1 





-2254.29 
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-337.24 
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-1954.76 


10 


-243.42 
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-1669.25 


11 


-169.45 
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-1407.21 


12 


-112.19 
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-1167.64 


13 


-68.88 
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-951.37 


14 


-37.48 
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-760.21 


15 


-16.55 
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-594.38 


16 


-4.86 


8 


-453.64 


17 


-0.50 



variable representation (DVR) [35] , and their binding en- 
ergies are presented in Table [I] The zero-point energy is 
found to be D /(hc) = 2254.29 cm- 1 . 



B. LiMg molecule 

To the best of our knowledge, there exist two experi- 
mental investigations of the LiMg ground state proper- 



ties [2"2l |2"3] . The study in Ref. [23 provides molecular 
constants that we can use for comparison with our calcu- 
lation. This work was based on resonant photo-ionization 
spectroscopy of LiMg in the gas phase and gives an ac- 
curate bond energy Do/(hc) = 1330 cm -1 , equilibrium 
separation R e = 5.9ao, and vibrational frequency uj e /{hc) 
= 190 cm -1 . A density functional study in Ref. [53] of 
LiMg gives D e /(hc) — 1371 cm -1 , R e = 6.04ao, and 
cu e /(hc) = 180 cm -1 . A MRCI approach with large con- 
tracted basis sets and multiple sets of polarization func- 
tions is used in Ref. [55J to give predictions for the LiMg 
ground state molecular constants D e /(hc) — 1611.28 
cm -1 , R e — 5.88ao, and uj e /(hc) = 183 cm -1 . The value 
for D e significantly differs from both Ref. [29] and exper- 
iment. 

Here we apply RCCSD(T) and UCCSD(T) meth- 
ods with two different basis sets to calculate the 
potential and dipole moment of the LiMg ground 
state. Augmented polarized core-valence (aug-cc- 
pCVQZ) basis sets with Li:(15s9p5d3flg)/[8s7p5d3flg] 
and Mg:(19sl5p6d4f2g)/[9s8p6d4f2g] functions [Ml [35] 
are used to calculate the ground state potential energy 
at the RCCSD(T) level of theory, where all electrons in 
both Li and Mg are correlated. Aug-cc-pv5z-DK basis 
sets with Li: (14s8p4d3f2glh)/[6s5p4d3f2glh] and Mg: 
(20sl4p4d3f2glh)/[7s6p4d3f2glh] [35] are applied in the 
UCCSD(T) calculation. The calculation is optimized for 
scalar relativistic calculations of the all-electron Douglas- 
Kroll Hamiltonian. A comparison with the molecular 
constants measured in Ref. [23] shows that our Do & n d 
R e values are in good agreement with this experiment. 
Figure [2] shows the independent potential and dipole mo- 
ment. 

In Table [LT] we present a complete list of N — bound 
states of the X 2 E + ground state potential of 7 Li 24 Mg ob- 
tained by a DVR calculation. The zero-point corrected 
energy is found to be Do/(/ic)=1330.05 cm -1 in agree- 
ment with experimental measurement |23j . 



TABLE II: All vibrational energies of the 7 Li 24 Mg X 2 E+ 
ground state potential for the N = rotational state. The 
potential is based on the RCCSD(T) calculation. 
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Energy/ (he) (cm 1 ) 
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Energy/(ft,c)(cm 1 





-1330.05 


10 


-216.21 
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-1164.35 


11 


-164.06 
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-1010.96 


12 


-120.16 
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-870.00 


13 


-83.91 
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-741.52 


14 


-54.87 
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-625.41 


15 


-32.62 
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-521.43 


16 


-16.78 
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-429.16 


17 


-6.89 
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-348.00 


18 


-1.75 
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-277.27 


19 


-0.12 
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FIG. 2: Top panel: electronic X 2 E + ground state potential 
of LiMg; bottom panel: permanent dipole moment, where a 
vertical dashed line indicates the equilibrium separation of 
the ground state potential; There the dipole moment is 0.44 
a.u., while the vibrationally averaged value for the v — level 
is 0.413 a.u.. 



C. LiCa molecule 

Currently, there is no reliable experimental data for 
the ground state potential of the LiCa molecule in the 
literature. However, preliminary experimental results on 
the electronic properties of the ground-state LiCa poten- 
tial based on high resolution spectroscopy were presented 
at the EGAS conference in 2010 [IT]. They gave the 
molecular constants D e /(hc) = 2607.8(100) cm -1 , R e — 
6.3415(5)a , and u> e /(hc) = 195.2 cnr 1 for 7 Li 40 Ca. Ear- 
lier theoretical predictions for the LiCa ground state give 
a dissociation energy of D e /(hc) = 2177.7 cm -1 and 
R e = 6.65ao using density functional theory [53]. Later 
a configuration-interaction calculation involving an ef- 
fective pseudo-potential [38] predicted D e /(hc) = 2355 
cm -1 and R e = 6.23ao- An all-electron ab initio study 
of LiCa in Ref. [15] reported D e /(hc) = 1935.7 cm Jl 
and R e — 6.44ao- Experimental two-photon ionization 
spectroscopy reported in the same work found that R e 
equals to 6.3865a . The substantial discrepancy in the 
reported characteristics and lack of the experimental data 
on the electronic structure of this molecule has initiated 
our study. 

An accurate treatment of the LiCa molecule with 
presumably weak bonding requires a massive inclusion 
of correlation effects. To satisfy this requirement we 



apply two coupled-cluster approaches, UCCSD(T) 
and UCCSDT, that allow us to introduce higher- 
order correlation effects. To probe the influence of 
basis sets on the LiCa characteristics we used two 
different basis sets for each atom. These are all- 
electron basis sets with quadruple-zeta quality. Basis 
sets (def2-QZVPP) for Li:(15s6p2dlf)/[6s4p2dlf] 
and Ca:(24sl8p6d3f)/[lls6p4d3f] functions are 
taken from Ref. [IT] . Basis sets (aug-cc-pCVQZ) 
with Li:(15s9p5d3flg)/[8s7p5d3flg] [34] and 
Ca:(25sl9pl0d4f2g)/[10s9p7d4f2g] [37 functions in- 
clude larger correlation expansions. 

We find that the molecular characteristics depend not 
only basis set but also on the number of correlated elec- 
trons n e i. For example, for a 3-electron correlation one 
2s 1 valence electron of Li and two 4s 2 valence electrons 
of Ca are correlated, whereas for 13-electron correlation 
all electrons of Li and ten 3s 2 3p 6 4s 2 electrons of Ca are 
included in the active space. We conclude that using 
the UCCSD(T) method and the more advanced basis set 
(aug-cc-pCVQZ) leads to an increase of D e /(hc) and d e 
by w 5%. For the UCCSDT method we expect similar 
trends. We believe that the UCCSDT calculation with 
the basis set (aug-cc-pCVQZ) provides the most accu- 
rate electronic potential and dipole moment. These are 
shown in Fig. [3] Moreover, this calculation is in a good 
agreement with the experimental molecular constants re- 
ported by Ref. [TT] . 

The upper panel of Fig. [3] shows the potential energy 
curve of the X 2 S + state of LiCa as a function of inter- 
atomic separation. The permanent dipole moment of this 
state as a function of R is presented in the lower panel. 
A complete list of the N — vibrational energies of the 
ground state is presented in Table |HT[ 



TABLE III: All vibrational energies of the 7 Li 40 Ca X 2 E+ 
ground state potential for the rotational state N — 0. The 
potential is based on the UCCSD(T)/aug-cc-PCVQZ calcu- 



lation. 


V 


Energy /(he) (cm 1 ) 
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Energy/(/ic) (cm 1 ) 





-2428.02 


14 


-442.18 
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-2232.09 


15 


-366.32 


2 


-2042.75 


16 


-298.77 


3 


-1861.87 


17 


-239.22 


4 


-1688.77 


18 


-187.31 


5 


-1524.31 


19 


-142.70 


6 


-1368.30 


20 


-105.03 


7 


-1221.05 


21 


-73.95 


8 


-1082.74 


22 


-49.08 


9 


-953.39 


23 


-30.02 


10 


-833.10 


24 


-16.33 


11 


-721.92 


25 


-7.41 


12 


-619.79 


26 


-2.47 


13 


-526.60 


27 


-0.42 
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FIG. 3: Top panel: electronic ground state 2 S^ /2 potential of 
LiCa; bottom panel: permanent dipole moment of the ground 
state of LiCa. The vertical dashed line indicates the equilib- 
rium separation of the ground state potential. At this sepa- 
ration the dipole moment is 0.044 a.u. The vibrationally av- 
eraged dipole moment with the v — wave function is 0.437 
a.u. 



FIG. 4: Top panel: electronic ground-state 2 E+ potential 
of LiSr; bottom panel: permanent dipole moment of the 
ground state of LiSr. The potentials and dipole moments are 
obtained by two different methods: UCCSDT (solid curve) 
and UCCSD(T) (dashed curve). The vertical dashed line in 
the bottom panel indicates the equilibrium separation of the 
ground state potential. The vibrationally averaged dipole mo- 
ment with the v = wave function is 0.111 a.u.; 



D. LiSr molecule 



To our knowledge there exist no data on the LiSr 
ground-state potential based on experimental observa- 
tion. We are only able to compare our results to the the- 
oretical results of Ref. [H] . The calculation of Ref. [5] uses 
an effective core and core polarization potentials that 
are incorporated in a configuration-interaction approach. 
Their molecular constants are R e — 6.57ao, D e /(hc) — 
2587 cm" 1 , u e /(hc) = 184.9 cm" 1 , B e /(hc) = 0.21 cm" 1 , 
and d e = 0.13 a.u. In the present study we use different 
computational methods to predict these characteristics. 

We perform calculations for LiSr that are similar to 
those for LiCa. We use UCCSD(T) and UCCSDT meth- 
ods with three different basis sets. The inner electrons of 
the closed ls 2 2s 2 ...3d 10 shells of Sr are chemically inac- 
tive and excluded from the correlation calculation. The 
ls 2 2s 1 electrons of Li and the ten 5s 2 4s 2 4p 6 electrons of 
Sr, are included in the all-order correlation calculation. 

In the first calculation we use the basis set 
def2-QZVPP with Li:(15s6p2dlf)/[6s4p2dlf] and 
Sr:(8s8p5d3f)/[7s5p4d3f] functions [41] and the 
Stuttgart ECP28MDF relativistic effective core po- 
tential (ECP) 40 . In a second calculation we use 
the basis set "[32]" with Li: (14s8p4d2f)/[7s6p4d2f] 



and Sr: (14sllp5d4flg)/[ 8s8p5d4flg] functions. The 
effective core potential is from Ref. [32] ■ Finally, 
we performed a calculation with a more-advanced 
all-electron basis set (aug-cc-pCV5Z) described 
by Li: (18sl2p7d5f3glh)/[10s9p7d5f3glh] and Sr: 
(23sl9pl2d4f2g)/[llsl0p7d4f2g] functions from [35]. 
The UCCSDT calculation is time consuming and has 
been performed for basis set "[32]" only. 

Results of the most advanced calculation by the 
UCCSDT method with the basis set "[32]" and the 
UCCSD(T) method with basis set (aug-cc-pCV5Z) are 
shown in Fig. [4] We see that the potential energy of 
the two calculations nearly coincide, whereas the dipole 
moment curves in the bottom panel are slightly shifted 
relative to each other. We believe that this difference is 
within the accuracy of our calculation. The largest dis- 
crepancy of about 7% with the results of Ref. [5] occurs 
at large internuclear separation, whereas at R e the dipole 
moment d e differs by 2% to 3% only. 

The DVR calculation of the bound state energy gives 
Do / (he) = 2275.59 cm" 1 and u/(hc) =182.2 cm" 1 . Ta- 
ble IV lists all N — bound-state energies for 7 Li 88 Sr. 
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TABLE IV: All vibrational energies of the 7 Li 88 Sr X 2 E+ 
ground state potential for the rotational state N = 0. The 
potential is based on the UCCSD(T) calculation with basis 
set [35] 



1000 
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Energy/ (he) (cm 1 ) 
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Energy/(/ic) (cm 





-2275.59 


15 


-397.78 
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-2099.56 


16 


-331.60 
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-1931.06 


17 


-272.33 
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-1769.35 


18 


-219.76 
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-1614.41 


19 


-173.66 
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-1466.66 


20 


-133.76 
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-1326.14 


21 


-99.83 
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-1192.91 


22 


-71.56 


8 


-1067.12 


23 


-48.65 
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-948.89 


24 


-30.80 


10 


-838.24 


25 


-17.64 


11 


-735.17 


26 


-8.73 


12 


-639.65 


27 


-3.42 


13 


-551.65 


28 


-0.86 
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-471.08 


29 


-0.04 
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E. LiYb molecule 

Recently, two extensive quantum-chemical calculations 
of the ground and excited states of the LiYb molecule 
have been reported by Zhang et al. 8J and by Gopaku- 
mar et al. [S] . The first calculation used both MRCI and 
UCCSD(T) methods. Scalar relativistic effects were in- 
cluded by the Douglas-Kroll Hamiltonian and an effective 
core potential. Their molecular constants for the X 2 £ + 
ground state of 7 Li 172 Yb are the following: R e = 6.681ao, 
D e /(hc) = 1577 cm" 1 , w e /(hc) = 147.36 cm" 1 . Different 
treatments of the Yb core electrons lead to different val- 
ues of the dipole moment. A RECP-pseudo-potential cal- 
culation gave d e = 0.022 a.u., while an all-electron DKH3 
calculation gave a smaller d e of 0.011 a.u. Results of a 
calculation by [5] were based on both the CASPT2 per- 
turbation theory and the CCSD(T) method. Relativis- 
tic effects were taken into account through the Douglas- 
Kroll-Hess Hamiltonian. Their best values of molecular 
constants are: R e = 6.669ao, D e /(hc) — 1421.96 cm -1 , 
u) e /(hc) = 135.54 cm" 1 . 

Our goal here is to perform an independent calcu- 
lation of the ground state potential and permanent 
dipole moment of LiYb using the UCCSD(T) method 
with another effective ECP28MWB pseudo-potential and 
with basis sets for Li: (15s6p2dlf)/[6s4p2dlf] and Yb: 
(14sl3pl0d8f6g)/[10s8p5d4f3g] functions from 03]. 

Results of our calculation are shown in Fig. [5] The 
upper panel of Fig. [5] shows the ground state potential 
of LiYb, which is significantly deeper than the potentials 
obtained in Refs. JSUH]- The permanent dipole moment 
is presented in the bottom panel. It is small and agrees 



FIG. 5: Top panel: electronic ground-state 2 E + potential of 
LiYb; bottom panel: permanent dipole moment of the ground 
state. The vertical dashed line in the bottom panel indicates 
the equilibrium separation of the ground state potential; 



well with the DKH3 value of Ref. [8]. 

Table W\ shows all N = bound-state energies for 
7 Li 172 Yb. 



II. SUMMARY 

The objective of our research has been a systematic 
investigation of the X 2 E + ground state properties of Li- 
group II molecules. To achieve this goal we applied 
powerful quantum chemistry UCCSD(T) and UCCSDT 
methods with large basis sets to calculate potential ener- 
gies and permanent dipole moments. Table |VI| displays 
molecular constants obtained in our calculations. 

Figures [6] and [7] and Table VI summarize our results. 
Figure [6] reveals a most striking irregularity in the disso- 
ciation energy and R e among the Li-group II molecules. 
The deepest and most shallow potentials belong to the 
LiCa and LiMg molecules, respectively. 

Figure [7] shows dipole moment functions for all 
molecules considered in our study. There is some con- 
sistency in the maximum value of d with the reduced 
mass of the molecule: heavier molecules have smaller 
dipole moment. The dipole moment averaged over the 
ground state vibrational wave function is always a few 
tenths of atomic units smaller than d e . The LiSr and 
LiYb molecule have the smallest dipole moments among 
all molecules in this study. 
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TABLE V: All vibrational energies of the 7 Li 172 Yb X 2 E+ 
ground state potential for the rotational state N — 0. The 
potential is based on the UCCSD(T) calculation with basis 

set m 



V 


Energy/ (he) (cm ) 


u 


Energy/ (he) (cm ) 





-2277. .12 


16 


-349.63 


1 


-2103.98 


17 


-289.41 


2 


-1938.15 


18 


-235.73 


3 


-1778.86 


19 


-188.38 


4 


-1626.11 


20 


-147.10 


5 


-1480.30 


21 


-111.67 


6 


-1341.45 


22 


-81.82 


7 


-1209.63 


23 


-57.26 


8 


-1084.99 


24 


-37.70 


9 


-967.64 


25 


-22.84 


10 


-857.62 


26 


-12.36 


11 


-754.90 


27 


-5.66 


12 


-659.49 


28 


-1.94 


13 


-571.36 


29 


-0.36 


14 


-490.44 


30 


-0.005 


15 


-416.57 







FIG. 6: Electronic ground state 2 E + potentials of Li-group II 
molecules and LiYb as a function of interatomic separation. 
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FIG. 7: Permanent dipole moment of the ground state of the 
Li-group II molecules and LiYb as a function of interatomic 
separation. Filled circes indicate dipole moment values at 
equiliblium separations of corresponding molecules. 



Table |VI| presents molecular constants calculated by 
the RCCSD(T), UCCSD(T) and UCCSDT methods us- 
ing several large basis sets available in the literature. 
Comparison of our most accurate (either UCCSD(T) or 
UCCSDT) calculations with existing experimental mea- 
surements of the molecular constants shows that we are 
able to reproduce the bond length (R e ) , harmonic fre- 
quency (w e ), and dissociation energy (D e ) within 5%, 
8%, and 6%, respectively. 

The dipole moment as a function R is calculated in 
terms of a charge transfer and induced charge transfer. 
These two contributions have opposite signs over a large 
range of internuclear separations. We conclude from our 
calculations that charge transfer from the Li-atom to 
the alkaline-earth atom prevails at all separations. This 
trend is reversed for LiYb at R > 7ao. 

We also highlight the importance of the correlation ef- 
fects in calculations of the alkali-alkaline-earth molecules. 
We have gathered numerical evidence to show that molec- 
ular constants of these molecules, such as dissociation en- 
ergy and permanent dipole moment of the ground state, 
are extremely sensitive to high order correlations. We 
believe that the best values of these characteristics are 
obtained by using the UCCSDT method, where correla- 
tion effects are taken into account at all orders. 

Finally, we want to emphasize that LiBe and LiMg 
are not likely candidates for ultracold studies because of 
the difficulty of cooling Be or Mg. However, LiCa looks 
quite promising, since a Ca BEC [45] has been made, 
and its dipole moment is comparable to the known RbCs 
case [35]. The 2 S + ground state gives an extra spin de- 
gree of freedom relative to the 1 £ + RbCs species, thus 



TABLE VI: Molecular constants for the X 2 E+ state of LiBe, 
LiMg, LiCa, LiSr, and LiYb obtained by the RCCSD(T), 
UCCSD(T) and UCCSDT methods. Here, 1 a =0.0529177 
nm is the Bohr radius and a dipole moment of 1 a.u. cor- 
responds to 1 eao, where e is the electron charge. Also h is 
Planck's constant and c the speed of light. 



Description 


Re 


W e / (Ac) 


D e /{hc) 


d t 




(o ) 


(cm" 1 ) 


(cm" 1 ) 


(a.u.) 


7 Li 9 Be 


UCCSD(T) 










aug-cc-pV5Z-DK 


4.873 


299.5 


2406 


1.41 


Experiment [32] 


4.894 








7 Li 24 Mg 


RCCSD(T) 










aug-cc-pCVQZ 


5.86 


174.4 


1417 


0.32 


UCCSD(T) 










aug-cc-pV5Z-DK 


5.87 


206.1 


1432 


0.44 


Experiment 23 


5.9 


190 


1423° 




7 Li 40 Ca 


UCCSD(T) 










def2-QZVPP 


6.410 


196.0 


2320 


0.438 


aug-cc-pCVQZ 


6.358 


205.5 


2460 


0.456 


UCCSDT 










aug-cc-pCVQZ 


6.357 


207.1 


2607 


0.440 


Experiment [TT] 


6.3415(5) 195.2 


2607.8(100) 


7 Li 88 Sr 


UCCSD(T) 










def2-QZVPP 


6.766 


167.0 


2165 


0.109 


basis of [42] 


6.712 


183.0 


2302 


0.117 


aug-cc-pCV5Z 


6.700 


182.2 


2367 


0.112 


UCCSDT 










basis of [42] 


6.711 


184.2 


2401 


0.096 


7 Li 172 Yb 


UCCSD(T) 










basis of [H] 


6.710 


181.5 


2289 


0.011 



more quantum control possibilities. While ultracold Li 
+ Sr and Li + Yb mixtures are possible for making 
the molecules, the dipole moments of the ground state 
molecules are rather modest. It remains to be deter- 
mined how best to make the molecules from the trapped 
atoms. One possibility is to make "preformed pairs" in 
a dual species optical lattice [UJ US] , and then use inter- 
combination line photoassociation followed by STIRAP 
[4"5] . Although all of these molecular species are expected 
to be highly chemical reactive in their ground states since 
the formation of the dimers is exotermic LiCa might have 
a large enough dipole moment to allow it to be protected 
against destructive collisions by aligning it with an elec- 
tric field in a ID optical lattice [U [50j [51] . 
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